Journal  of  Power  Sources  196  (2011)  8525-8530 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


in 

SllbaidtS 


Bi-modal  water  transport  behavior  across  a  simple  Nation  membrane 

Ziheng  Zhang31,  Keith  Promislow5,  Jonathan  Martinc,  Haijiang  Wangc,  Bruce  J.  Balcom3’* 

a  MRI  Centre,  Department  of  Physics,  University  of  New  Brunswick,  Fredericton,  New  Brunswick,  Canada  E3B  5A3 
b  Department  of  Mathematics,  Michigan  State  University,  East  Lansing,  MI  48824,  USA 

c  Institute  for  Fuel  Cell  Innovation  (IFCI),  National  Research  Council  Canada  (NRC),  4250  Wesbrook  Mall,  Vancouver,  BC,  Canada  V6T  IW5 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  9  June  2011 
Accepted  10  June  2011 
Available  online  22  June  201 1 


Keywords: 

PEM  fuel  cell 
Water  management 
Nation 

Mathematical  model 

Nuclear  magnetic  resonance  imaging 

MRI 


The  development  of  predictive  mathematical  models  for  water  management  in  polymer  electrolyte  mem¬ 
brane  fuel  cells  requires  detailed  understanding  of  water  distribution  and  water  transport  across  the 
Nation  layer.  The  anisotropic  microstructure  of  Nation  suggests  the  measurement  of  water  content  and 
mass  transport  should  be  along  the  fuel  cell  functional  direction,  i.e.  across  the  membrane.  Non-invasive, 
high  resolution,  microscopy  measurements  of  this  type  are  very  challenging. 

We  report  here  the  calibration  of  a  minimal  mathematical  model  for  diffusive  water  transport  in  Nation 
against  data  from  high-resolution  water  content  maps  determined  with  a  new  magnetic  resonance  imag¬ 
ing  methodology  developed  for  this  purpose.  A  mock  fuel  cell  was  designed  to  permit  well-controlled 
wetting  and  drying  boundary  conditions.  With  no  chemical  potential  driving  force  involved,  we  assume 
the  water  transport  behavior  will  be  dominated  by  diffusion.  Moreover  we  show  that,  in  this  context, 
our  model  is  mathematically  equivalent  to  the  traditional  permeation  models  based  upon  saturation 
dependent  pressure  gradients  via  a  capillary  pressure  ansatz. 

The  non-linear  equilibrium  water  distribution  across  the  Nation  membrane  measured  in  this  work 
suggests  a  bi-modal  diffusivity.  The  model  constructed  associates  distinct  transport  behaviors  to  water 
contents  above  and  below  a  critical  threshold,  consistent  with  a  rearrangement  of  a  micro-structural 
pore  network.  The  experimental  observation  and  the  model  prediction  agree  with  the  primary  features 
of  Weber’s  model  of  Nation,  which  predicts  distinct  modes  of  transport  for  hydration  fronts  traversing 
the  through-plane  direction  of  the  membrane. 

©  201 1  Published  by  Elsevier  B.V. 


1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cells,  as  power 
conversion  engines,  generate  useful  voltage  from  electrochem¬ 
ical  energy  for  automotive  and  stationary  power  applications 
with  high  efficiency  and  benign  emissions  [1-3].  The  heart  of 
PEM  fuel  cells  is  the  organic  polymer  electrolyte  membrane, 
e.g.  Nation,  functioning  simultaneously  as  an  electrolyte  conduc¬ 
tor  and  a  gas  separator.  The  efficiency  of  polymer  electrolyte 
membranes  for  fuel  cell  applications  is  dictated  by  the  proton 
conductivity,  which  is  determined  by  the  proton  mobility  and 
proton  concentration  in  the  aqueous  phase  of  the  membrane 
[11- 

Water  management,  maintaining  an  appropriate  water  balance 
within  the  membrane  through  control  of  the  relative  humidity  (RH) 
of  the  input  reactant  gases,  is  critical  to  ensuring  optimal  operation 
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of  the  PEM  fuel  cell  [4,5].  Excess  water  produced  in  PEM  fuel  cells  at 
the  cathode  can  limit  reactant  gas  transport.  Too  little  water  at  the 
anode  can  cause  the  membrane  to  dry  out  and  lose  its  conductivity. 
Reliable  water  management  requires  a  detailed  understanding  of 
the  transport  mechanisms  of  water  through  the  Nation,  which  are 
still  not  well  established. 

Many  models  propose  diffusion  as  the  basic  transport  mech¬ 
anism  for  water  in  Nation  [4],  while  others  suggest  hydraulic 
permeation  [6-10].  Many  current  modeling  efforts  suppose  a 
single-mode  transport  mechanism  for  water  through  the  Nation, 
and  hence  predict  a  smooth  water  distribution  across  the  mem¬ 
brane.  However,  recent  experimental  evidence  shows  that  a  dry 
anode  can  be  maintained  while  the  wetter  cathode  has  liquid  water 
present  during  the  fuel  cell  operation  [1,9]. 

Moreover,  Weber  et  al.  [11,12]  have  proposed  a  model  that  pro¬ 
poses  various  transport  modes  based  upon  pore-expansion  due 
to  liquid  pressure  in  connective  nanopores  within  the  ionomer. 
Depending  upon  the  saturation  level,  their  model  has  co-existing 
diffusive  and  permeation-based  transport.  In  contrast  to  other  the¬ 
ories,  this  model  incorporates  a  degenerate  transport  mechanism, 
which  indicates  distinct  transitions  in  water  profiles  without  large 
fluxes  [11,12]. 
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In  an  operating  fuel  cell,  the  proton  conductivity,  water  dif¬ 
fusion  coefficient,  and  electro-osmotic  drag  coefficient  all  vary 
as  functions  of  water  content  in  the  membrane  [13].  Direct 
measurements  of  water  content  in  Nation  are  hence  of  great 
importance  in  developing  appropriate  mathematical  models,  and 
many  related  studies  have  been  reported  in  the  last  decade  [14]. 
However,  the  anisotropic  microstructure  of  Nation  [15-17]  sug¬ 
gests  the  importance  of  measurements  of  water  content  along  the 
cell’s  functional  direction,  i.e.  across  the  membrane,  which  is  dis¬ 
tinct  from  information  gathered  from  lengthwise  measurements 
[18]. 

At  present,  optical  visualization  [19-23],  neutron  beam  absorp¬ 
tion  [24-29],  and  magnetic  resonance  imaging  [13,30-34]  are  three 
major  non-invasive  techniques  utilized  to  monitor  the  water  con¬ 
tent  variation  inside  PEM  fuel  cells.  Optical  methods  depend  on 
the  transparency  of  the  detecting  materials.  However,  Nation  is 
semi-transparent,  and  will  not  permit  quantitative  water  measure¬ 
ments  as  a  function  of  depth  across  the  membrane.  Neutron  beam 
methods  feature  a  high  sensitivity  to  water  content;  however  they 
suffer  from  poor  spatial  resolution  and  difficulty  with  coordinate 
registration. 

Magnetic  resonance  imaging  is  a  promising  technique  with  the 
potential  for  high  resolution  measurement  [35-39].  From  the  lit¬ 
erature,  the  majority  of  MRI  approaches  to  this  problem  are  based 
on  frequency  encoding  spin  echo  techniques,  which  are  highly  sen¬ 
sitive  to  the  local  magnetic  susceptibility  variation  and  often  yield 
images  with  artifacts.  In  general,  removing  these  image  artifacts 
requires  a  wide  bandwidth  signal  which  leads  to  a  decreased  signal 
to  noise  ratio  (SNR)  [33].  Moreover,  the  inherent  signal  contrast, 
which  depends  on  the  spin-spin  relaxation  time,  T2 »  renders  the 
frequency  encoded  images  non-quantitative,  unless  T2  mapping  is 
undertaken  [38,39]. 

In  this  study,  we  report  the  measurement  of  the  water  content 
distribution  across  a  Nation  1110  membrane  in  a  mock  fuel  cell, 
purposely  designed  to  study  water  content  variation  across  a  Nation 
membrane  with  one  side  equilibrated  by  liquid  water  and  the  other 
exposed  to  water  vapor  at  well  controlled  RH.  In  the  absence  of  the 
chemical  reactions  present  in  an  operational  fuel  cell,  our  experi¬ 
mental  set-up  eliminates  the  electromotive  driving  force  for  proton 
transfer,  and  focuses  on  water  transport  issues.  The  results  reveal  a 
surprising  family  of  quasi-equilibria  of  water  distributions.  In  each 
time  resolved  image,  the  water  content  profile  through  the  water 
equilibrated  side  of  the  Nation  is  relatively  flat,  with  a  slight  down¬ 
ward  tilt  toward  the  vapor  equilibrated  side,  however  at  a  critical 
value  of  the  water  content  the  slope  of  the  water  profile  increases 
significantly.  We  attribute  this  non-uniform  behavior  to  a  phase- 
transition  in  the  microstructure  of  the  underlying  water  network 
[40], 

This  observation  is  similar  to  the  observation  of  wet- 
cathode/dry-anode  steady  state  in  an  operating  fuel  cell  reported 
in  several  previous  studies,  which  cannot  be  well  described 
by  non-degenerate  transport  models  [41,42].  In  this  study, 
motivated  by  a  reduction  of  the  Weber-Newman  model,  we 
propose  a  minimal  model  for  water  transport,  using  a  self- 
diffusivity  for  water  that  depends  sensitively  upon  water 
concentration  at  a  particular  threshold.  This  model  describes 
steady  state  behavior,  with  slow-relaxation  or  aging  effects 
ignored. 

The  quantitatively  excellent  agreement  between  the  experi¬ 
mental  measurements  and  the  model  results  support  the  idea 
of  bifurcation  in  water  network  structure,  as  parameterized  by 
average  water  content,  leading  to  a  bi-modal  water  transport  mech¬ 
anism.  Moreover  the  transition  threshold  of  our  calibrated  model  is 
in  excellent  agreement  of  the  liquid-vapor/vapor  phase  transition 
predicted  by  the  Weber-Newman  model. 


Fig.  1.  Schematic  of  the  mock  fuel  cell  with  wetting  and  drying  boundary  conditions. 
The  Nafion  layer  (b)  is  held  between  two  perforated  plates  (a,  c),  whose  structure  is 
illustrated  in  the  inset  at  right.  Plate  (a)  is  fabricated  from  copper  plated  PC  board. 
Plate  (c)  is  made  of  G10  fiberglass.  A  water  reservoir  sits  on  top  of  plate  (a)  with  a 
drying  chamber  underneath  plate  (c),  also  made  from  G10  fiberglass.  Beneath  plate 
(c),  a  circular  RF  coil  with  1.8  cm  diameter  (d)  was  installed.  The  entire  apparatus 
was  held  together  by  4  brass  screws  with  a  controlled  torque  applied. 


2.  Experimental 

2A.  Mock  fuel  cell 

As  shown  in  Fig.  1 ,  the  mock  fuel  cell  chamber,  3  cm  x  3  cm,  was 
shaped  from  G10  fiberglass  (McMaster-Carr,  Cleveland,  OH).  The 
water  reservoir  on  top  of  the  Nafion  layer  saturates  the  sample, 
when  desired,  to  ensure  a  fully  wet  boundary  condition.  On  the 
other  side,  the  gas  chamber  underneath  maintains  a  constant  dry¬ 
ing  boundary  condition  through  a  dry  air  flow,  with  the  RH  12%. 
A  perforated  PC  board  between  the  reservoir  and  membrane  was 
employed  to  remove  the  signal  from  the  water  reservoir  through 
RF  screening  [38]. 

All  perforated  plates  were  1.65  mm  in  thickness.  The  water 
reservoir  and  gas  chamber  frame  were  6.35  mm  thick.  The  assembly 
was  embedded  in  a  1  cm  thick  Teflon  sheet  and  placed  in  a  copper 
covered  Plexiglas  tube. 

2.2.  Sample  preparation 

The  sample  film,  Nafion  1110,  250  p,m  thick,  was  purchased 
from  Ion  Power  (New  Castle,  DE).  The  Nafion  1110  was  cut  to  a 
2  cm  diameter  circle,  and  pretreated  by  heating  in  3-5wt.%  H202 
solution  at  80  °C  for  1  h.  The  sample  was  then  rinsed  in  deionized 
water  at  80  °C  for  1  h,  then,  heated  in  0.5  M  H2S04  at  80  °C  for  1  h. 
The  sample  has  finally  cleaned  by  prolonged  heating  in  deionized 
water  at  80  °C  for  1  h,  to  remove  any  possible  chemical  residues. 
The  basic  procedure  was  based  on  the  work  of  Jalani  [43  ]  and  Moore 
[44]. 

2.3.  Magnetic  resonance  imaging 

The 1 H  magnetic  resonance  imaging  measurements  were  imple¬ 
mented  on  a  2.4 T  (100  MHz)  horizontal  bore  superconducting 
magnet,  with  a  Tecmag  Apollo  as  the  imaging  console  [38].  The 
RF  probe  was  a  custom  built  single  channel  surface  coil.  The  RF  res¬ 
onator  was  formed  from  copper  tape,  cut  into  a  circular  shape,  2  cm 
in  diameter,  affixed  to  the  surface  of  the  gas  chamber,  shown  in 
Fig.  1.  The  RF  resonator  was  tuned  to  100  MHz  with  two  capacitors 
in  a  standard  series-parallel  tank  circuit. 

A  specially  developed  MRI  technique  for  thin  film  depth  imag¬ 
ing,  Double  Half  K-space  Spin  Echo  Single  Point  Imaging  (DHK  SE 
SPI)  [37],  was  employed  to  obtain  ID  water  content  profiles.  As 
a  pure  phase  encoding  measurement,  it  minimizes  the  potential 
image  distortion  from  conductive  materials  in  a  fuel  cell  [45,46]. 
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From  previous  studies  [38,39],  it  is  well  known  that  the  water  T2  in 
Nation  varies  with  the  water  content  and  pore  level  microstructure. 

The  MRI  measurement  averages  the  cross  section  of  the  sam¬ 
ple  while  encoding  along  the  depth  direction  y  [38].  The  spatially 
resolved  image  intensity  is  simply, 

s(y)  =  p0(y)  •  exp  (-^))  (i) 

where  n  is  the  echo  number,  and  TE  is  the  echo  time.  S(y)  is  the 
image  intensity  defined  by  the  product  of  the  water  content,  p0(y)* 
and  T2  weighting  in  the  exponential  term. 

The  MRI  measurement  simultaneously  maps  T2,  permitting  one 
to  directly  determine  water  content  variation.  By  fitting  each  pixel 
in  the  profile  series  to  Eq.  (1 ),  it  is  possible  to  separate  the  spatially 
resolved  T2  and  the  local  density,  po{y),  corresponding  to  the 
water  content.  From  previous  studies  [38,39],  the  water  T2  in  Nation 
exhibits  a  single  exponential  decay,  which  is  assumed  in  the  T2 
fitting  for  all  measurements  in  the  current  study. 

The  nominal  resolution  of  the  images,  and  the  resultant  water 
content  profiles,  was  7.4  fxm  with  an  image  field  of  view  (FOV) 
of  490  p,m.  A  series  of  ID  profile  were  continuously  acquired  to 
temporally  resolve  the  experiment,  with  a  frame  rate  of  one  image 
every  610  s. 

3.  Results  and  discussions 

3.1.  Wetting  and  drying  measurements 

Unlike  operating  fuel  cells,  the  mock  fuel  cell  does  not  form  a 
complete  electrical  circuit;  it  has  no  catalyst,  no  gas  diffusion  lay¬ 
ers  and  no  chemical  reactions.  The  mock  fuel  cell  does  however 
provide  well-controlled  boundary  conditions.  The  evolution  of  the 
water  content  across  the  Nation  membrane  was  investigated  with 
a  wetting/drying  boundary  condition. 

As  illustrated  in  Fig.  1 ,  a  water  reservoir  on  the  top  of  the  Nation 
layer  was  employed  to  maintain  saturation  and  a  dry  gas  flow 
underneath  was  employed  for  drying.  In  this  case,  we  anticipate 
water  flow  into  the  top  of  the  membrane  to  reach  a  steady  state 
with  drying  from  below.  The  wetting/drying  experiment  mimics 
possible  operating  conditions  of  a  dry-gas-feed  fuel  cell  in  which 
the  cathode  is  exposed  to  liquid  water  while  the  anode  is  relatively 
dry. 

The  experiment  duration  was  approximately  20  h  with  more 
than  100  multi-echo  profiles  acquired,  at  approximately  10  min 
intervals  [38].  The  steady  state  lasted  for  approximately  1 8.5  h  until 
the  water  reservoir  emptied. 

Fig.  2  reports  4  averaged  water  content  profiles  of  this  wet¬ 
ting/drying  experiment.  Each  is  averaged  from  the  isolated  water 
content  profile  of  20  measurements,  which  overall  span  3.4  h.  The 
water  reservoir  is  at  the  left  in  the  profiles  with  the  drying  chamber 
at  right.  It  is  hard  to  determine  p0(y)  at  the  profile  edges,  because 
the  second  or  subsequent  echo  signals  of  these  points  are  usu¬ 
ally  too  small,  which  substantially  complicates  the  data  fitting.  The 
Nation  sample  is  located  from  140  pan  to  360  p,m,  with  the  thick¬ 
ness  approximately  220  |jim. 

The  images  clearly  show  a  systematic  variation  in  water  content 
across  the  membrane.  High  water  content  shown  on  the  wetting 
side  (left)  indicates  a  well  maintained  water  content  close  to  the 
wetting  boundary.  Moving  from  left  to  right,  toward  the  drying 
side  of  the  water  content  profile  in  Fig.  2,  one  observes  an  ini¬ 
tially  shallow  decrease  of  water  content.  At  a  critical  point  in  all 
the  water  content  profiles,  we  observe  a  sharp  decrease  in  water 
content  moving  toward  the  drying  side. 

Each  of  the  four  water  content  profiles  of  Fig.  2  is  consistent  with 
the  bi-modal  transport  model,  however  they  also  exhibit  an  “aging” 
effect  in  which  the  water  content  decreases  over  many  hours.  This 


Fig.  2.  4  averaged  water  content  profiles.  Experiments  (•),  (O).  (▼),  (□)  corre¬ 
spond  to  the  average  from  the  experimental  time  range  of  0-3.4, 3.4-6.8,  6.8-10.2, 
10.2-13.6  h.  Water  content  uniformly  decreased  over  20  h,  with  a  clear  drying  front 
observed  at  the  drying  side  (right).  The  typical  experimental  uncertainty  due  to 
fitting  is  illustrated  by  the  representative  error  bar  in  the  top  profile  at  the  far  left. 


aging  effect  is  quantified  via  systematic  changes  in  the  spatial  res¬ 
olution  of  the  T2  lifetime  [38]  and  suggests  a  slow  evolution  of 
the  microscopic  Nation  port  structure  with  a  time  constant  of  tens 
of  hours,  as  suggested  in  the  literature  [47-51].  This  aging  effect 
implies  that  the  threshold  Xt  for  the  bi-modal  transport  model  may 
depend  upon  quantities  beyond  the  scope  of  our  model,  including 
the  past  history  of  the  Nafion  water  concentration. 


3.2.  Numerical  model 

A  medium  with  a  uniform,  isotropic  microstructure,  and  a 
constant  diffusivity,  would  predict  a  linear  variation  of  water  con¬ 
centration  between  the  two  boundaries.  This  is  in  fact  a  classic 
hallmark  of  diffusion  [52].  However,  the  nonlinear  water  distribu¬ 
tion  observed  in  our  measurement  suggests  that  a  mechanism  more 
complicated  than  simple  diffusion  is  controlling  the  mass  transport 
across  the  Nafion  membrane. 


3.2  A.  Reduction  of  the  Weber-Newman  model 

A  key  aspect  of  the  wetting/drying  data  in  Fig.  2  is  the  bi- 
modal  nature  of  the  water  transport.  Previous  models  for  water 
transport  in  Nafion  have  proposed  mechanisms  for  bi-modal  and 
even  tri-modal  mechanisms.  We  focus  on  the  model  presented  in 
[11,12]  which  combines  transport  models  for  vapor,  vapor-liquid, 
and  liquid  equilibrated  Nafion.  The  experimental  data  is  con¬ 
sistent  with  a  bi-modal  subset  of  this  model,  incorporating  the 
vapor  and  vapor-liquid  transport  mechanisms.  In  the  absence  of 
proton  transport,  the  analysis  below  shows  that  the  two-mode, 
Weber-Newman  model  is  equivalent  to  a  diffusive  transport  model 
in  which  diffusivity  depends,  sensitively,  upon  water  volume  frac¬ 
tion. 

The  Weber-Newman  model  is  presented  most  succinctly  in 
Eqs.  (1.1 )— (1.28)  of  [12],  which  we  summarize  below.  Taking  the 
transport  to  be  entirely  in  the  through-plane  direction,  x,  and  at 
equilibria,  the  vapor-liquid  model  prescribes  the  spatially  constant 
ionic  current,  /,  as 

1  =  ~S  TX0  +  SPl)  “  (1  - ^  fa  ^ ^  &Pl)  (2) 
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where  0  is  the  electric  potential  and  pL  is  the  liquid  pressure,  and 
S,  the  liquid  saturation  level,  is  an  affine  function  of  water  content, 
kt 


X  —  Xy 

XI  —  Xy 


0) 


In  particular  S  increases  from  0  to  1  as  the  water  content  varies 
from  its  saturated  vapor  level  X  =  XV  =  XV{T)  to  the  liquid  level 
X  =  Xl  =  22  waters/acid  group.  The  water  flux  N  is  also  spatially  con¬ 
stant  and  is  prescribed  by 


N -s  ('>■'  s4 + "“aM <’  - 51  (a-  s' »■ +fes»)  • 

(4) 

where  the  dependence  of  the  transport  coefficients  cr^X ,  T) 

and  pij(X,  T )  upon  water  content  and  temperature  can  be 
inferred  from  [12].  Moreover  the  liquid  pressure  is  related  to  the 
saturation  level  via  an  empirical  relation  based  upon  assump¬ 
tions  of  hydrophobicity  (cos(0)<O)  and  a  particular  pore  size 
distribution, 


Fig.  3.  “Bi-modal”  diffusivity  proposed  as  a  function  of  water/acid  group  number. 
The  threshold  transition  point  is  represented  by  Xt,  with  the  transition  region  rep¬ 
resented  by  e. 


S 


1 

2 


^1  -  erf 


/ln[-(5/2)ycos(e)]-lnpA\ 
V  0.3V2  ) ) 


(5) 


This  function  is  monotonic  in  pL  >  0,  and  hence  may  be  inverted 
to  slave  the  liquid  pressure  to  the  saturation  level.  However  the 
saturation  level  in  turn  is  a  function  of  water  density  X ,  which  yields 
a  capillary  pressure  relation, 


Pl=PlP), 


(6) 


for  water  fractions  X  e  [XV,XL],  in  the  vapor-liquid  regime. 

In  the  wetting/drying  experimental  configuration,  there 
is  no  ionic  current,  1  =  0  in  Eq.  (2),  so  that  the  electric 
potential  gradient  must  balance  against  the  water  pressure 
gradient, 


—  0-  5<7i2  +  (1  -  S)<722  d 

dx  ~  Sau  +  (1  -  S)<j22  dxPL 


(7) 


Combining  this  expression  with  the  dependence  of  the  liquid  pres¬ 
sure  upon  X  we  may  write  the  flux  as 

N  =  -Dvl(),,T )^X,  (8) 

where  the  effective  diffusivity  in  the  vapor-liquid  regime,  accord¬ 
ing  to  Eqs.  (4)  and  (6),  takes  the  form 


infinite.  This  fact  is  a  direct  consequence  of  the  choice  of  the  func¬ 
tional  relation,  Eq.  (5),  which  leads  to  the  corresponding  capillary 
pressure  relation,  Eq.  (6).  Indeed,  since  d/dpL  S  approaches  zero  as  S 
approaches  zero,  the  inverse  relation,  pL  =  Pl(S),  must  obtain  an  infi¬ 
nite  slope  at  S  =  0.  Equally,  the  derivative  of  the  capillary  pressure, 
which  is  a  factor  in  the  water  diffusivity,  satisfies  d/dXpL  (Ay)  =  oo. 
It  seems  more  plausible,  and  is  consistent  with  the  data  presented 
in  Fig.  2,  that  the  effective  water  diffusivity  should  increase  mono- 
tonically  with  increasing  water  content,  and  indeed  should  increase 
rapidly  with  a  transition  of  water  content  from  the  vapor  to  the 
vapor-liquid  regimes,  as  illustrated  in  Fig.  3.  To  this  end,  and  as 
is  consistent  with  the  level  of  resolution  afforded  by  the  data, 
we  retain  the  dominant  features  of  the  Weber-Newman  model, 
namely  a  bi-modal  transport  mechanism,  Fig.  3,  that  shifts  sensi¬ 
tively  at  the  transitional  water  content,  X  =  Xy,  but  simplify  the  flux 
relation  to 

N  =  ~D(X)^X,  (11) 

where  the  effective  diffusivity,  given  by  Eq.  (4),  is  strictly  increasing 
with  water  content,  and  exhibits  a  sharp  transition  from  a  linear 
dependence  upon  X  in  the  vapor  regime,  to  a  significantly  higher 
constant  value  in  the  vapor-liquid  regime  at  a  transition  value  Xt. 

3.2.2.  The  minimal  “ bi-modal  diffusivity"  model 

A  schematic  diagram  of  the  numerical  simulation  is  illustrated 
in  Fig.  4,  with  the  boundary  conditions  specified  on  the  left  and 


DvdK  T)  =  -  (^^(l-siS^n  +  0  “  S)ftl }  “  (S^12  +  0  “  S)fe)) 


(9) 


For  water  contents  X<XV,  the  liquid  saturation  S  is  considered 
to  be  zero,  and  the  transport  is  determined  from  the  gradient 
of  the  electric  potential  0  and  the  chemical  potential  /x0.  The 
zero-current  reduction  slaves  the  electric  potential  and  chem¬ 
ical  potential  gradients,  and  yields  a  simple  flux  relation  of 
form 

N  =  -«„(*,  =  -DV(K  T)^K  (10) 

where  the  coefficient  av  is  defined  in  Eq.  (1.15)  of  [12]  while  the 
chemical  potential  can  be  viewed  a  function  of  X  and  T  [  11],  which 
motivates  the  definition  of  the  effective  vapor-equilibrated  diffu¬ 
sivity,  Dy  =  av  d/dX  ii0. 

An  unintuitive  feature  of  the  Weber-Newman  model  is  that  as 
X  decreases  to  the  lower  limit,  Xv ,  of  the  vapor-liquid  regime,  the 
vapor-liquid  diffusivity  DVL  not  only  increases  but  even  becomes 


right  sides.  As  the  water  transport  problem  considered  in  this  study 
is  only  in  the  depth  direction,  possible  mass  transport  along  the 
membrane  surface  is  neglected.  The  water  concentration  is  denoted 


Water 


Gas  at 
RH=r0 


Y=L  Y=0 

Fig.  4.  A  schematic  diagram  of  the  proposed  mathematic  model.  The  thickness  of 
the  sample  is  I,  with  the  depth  direction  along  Y  direction,  according  to  the  axis 
system  in  Fig.  1.  The  drying  boundary  condition,  represented  by  r0,  is  listed  at  the 
right  side  of  the  sample  diagram,  with  the  wetting  boundary  condition,  saturated 
water,  listed  at  the  left  side. 
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by,  A  =  A(x),  in  dimensionless  units  of  water/acid  group.  The  actual 
water  concentration  is  a  •  A,  where  a  is  the  acid  density,  approxi¬ 
mately  1200  moles  m-3. 

Our  minimal  model  is  formulated  as 

'iH  I)"0  0<*<L 

<  D(\)^m=+y(m-Km .  (12) 

D(A)  •  =  -x(A(I)  -  A*(I)) 

<  ax 

where  y  is  the  mass  transport  coefficient;  D  is  the  water  diffusivity, 
which  depends  upon  water  content;  A *(RH)  represent  the  equilib¬ 
rium  membrane  water  content  as  function  of  external  RH,  which 
is  equivalent  to  the  water  vapor  activity.  Through  an  empirical 
formula,  A*(RH)  =  0.043  +  1 7.81  RH-  39.85  RH2  +36- RH3  [53,54], 
A *(RH)  is  calculated. 

In  our  proposed  “bi-modal”  diffusivity  model,  Fig.  3,  the  diffusiv¬ 
ity  increases  slowly  in  a  linear  manner  from  0  water/acid  number  to 
At,  but  jumps  from  a  concentration  near  At  to  a  substantially  higher 
value,  when  the  concentration  is  greater  than  At. 

This  diffusivity  coefficient  is  described  via  the  formulae, 

(  D  =  Dm  •  A  •  Hm(A)  +  Dp  •  HP(A) 

|  HP(X)  -  1  +  tanh^A  ~  1  ~  tanh((^  —  A.f)/g)  ’ 

(13) 


where  s  controls  the  width  of  the  diffusivity  transition  region,  while 
Dm  and  Dp  are  constants  calibrated  to  the  data. 

To  resolve  the  model,  Eq.  (2),  we  observe  that  the  flux 
N=-a  •  D(A)  •  dA/dx  is  constant  in  space  so  that  the  model  can  be 
rearranged  as, 


D(A). 


dA 

dx 


N 

a 


<  m  =  U0)-fL  . 

u  •  y 


A(L)  =  A*(L)  + 


N 

^7 


(14) 


We  introduce  the  cumulative  function  of  the  diffusivity,  D( A)  = 
fo  D( A)  •  dA,  so  that  d/dx(D( A))  =  -N/a,  and  hence, 

D(V  =  D0 (15) 


where  D0  =  D(A(0))  is  an  unknown  constant.  We  solve  for  A  in  terms 
of  N,  D0  as  X  =  Dt-DfDo  -  N  ■  x/a). 

To  determine  the  flux  N,  we  interpret  Eq.  (15)  from  the  bound¬ 
ary  conditions  and  obtain,  D(A(L))  -  D(A(0))  =  -N  •  L/a.  With  the 
combination  of  Eqs.  (12)-(14),  a  function  ofNis  obtained, 

D{K(L)  +  ^l)  ~  D  (K(0)  ~  Fl)  =“7  (16) 

One  is  therefore  able  to  solve  Eq.  (16)  for  N  through  New¬ 
ton’s  Method  and  determine  D0  from  the  relation,  D0  =  D(A(0))  = 
D(A^(0)  -  N/a  - 1). 


3.2.3.  Comparison  of  the  model  prediction  with  the  experiment 
We  attribute  the  overall  minor  decrease  of  the  water  content 
profiles  with  time,  Fig.  2,  as  being  due  to  the  slow  micro-structural 
rearrangement  of  Nation.  Since  this  is  a  minor  effect  with  a  con¬ 
sistent  shoulder  structure  among  the  four  profiles,  we  average  the 
four  water  content  profiles  to  remove  the  transient  aging  behav¬ 
ior  to  yield  steady-state  experimental  data,  which  can  be  used  for 
comparison  to  model  predictions,  as  illustrated  in  Fig.  5.  The  model 
curve,  Fig.  5,  is  computed  with  the  drying  and  wetting  boundary 


Fig.  5.  Comparison  of  the  computed  model  curve  (-)  against  the  wetting/drying 
experimental  data  (•),  an  average  of  the  four  profiles  shown  in  Fig.  2.  A  satisfying 
match  is  observed  which  demonstrates  the  correctness  of  the  proposition  of  a  bi- 
modal  diffusivity. 


conditions  set  as  0.12  and  1  according  to  the  RH  of  the  dry  air  flow 
and  the  water. 

The  difficulty  of  proper  T2  fitting  of  the  edges  of  the  T2-weighted 
profiles  yields  an  incomplete  water  content  profile.  As  a  result  the 
two  edges  of  the  profile  are  not  matched  to  the  boundary  condi¬ 
tions.  Hence,  in  this  study  the  experimental  water  profile  is  slightly 
rescaled  to  make  the  end  points  correspond  to  two  empirically 
estimated  RH  values.  Since  this  change  does  not  affect  the  pattern 
of  water  distribution  and  the  inherent  variation  of  diffusivity,  the 
compatibility  of  the  proposed  model  with  the  experimental  data  is 
therefore  maintained. 

Our  transition  value,  At  =  16  at  T=  20  °C,  is  in  excellent  agree¬ 
ment  with  the  Weber-Newman  value  of  Ay  for  an  unconstrained 
membrane,  as  extrapolated  from  the  data  presented  in  Fig.  3  of  Ref. 
[12].  Moreover,  the  satisfying  match  of  the  model  prediction  and 
experiment  observed  in  Fig.  5  supports  our  basic  assumption  of  a 
bi-modal  diffusivity.  In  addition,  this  match  indicates  the  existence 
of  a  critical  water  concentration  where  a  dramatic  change  in  diffu¬ 
sion  occurs,  which  also  suggests  a  threshold  may  exist  across  the 
Nation  in  an  operating  fuel  cell. 

To  our  knowledge,  this  experimental  data  afford  the  first  direct 
observation  of  such  a  nonlinear  water  distribution  through  a  Nation 
membrane  due  entirely  to  water  transport.  The  possible  underpin¬ 
ning  mechanisms  are  a  matter  of  speculation  at  this  time.  A  possible 
mechanism  could  include  a  phase  transition  in  the  water  pore  net¬ 
work  that  may  occur  at  a  critical  water  volume  fraction,  which 
in  turn  engenders  a  shift  in  either  diffusivity  or  permeability.  A 
full  characterization  of  the  nature  of  water  transport  and  water 
network  re-organization  is  the  subject  of  ongoing  work. 


4.  Conclusion 

This  paper  describes  measurement  of  the  water  content  vari¬ 
ation  across  a  Nation  membrane  inside  a  mock  fuel  cell,  under 
a  wetting/drying  boundary  condition.  The  water  content  profiles 
of  the  wetting/drying  measurement  present  a  behavior  that  is 
inconsistent  with  uniform  mass  transport  in  Nation.  A  minimal 
mathematic  model,  with  an  assumption  of  a  bi-modal  water  trans¬ 
port  mechanism  made  from  a  reduced  Weber-Newman  model,  is 
constructed  according  to  the  experimental  observations,  with  the 
subsequent  model  predictions  providing  excellent  match  to  the 
data. 
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Measurements  of  this  type  show  the  great  potential  of  ex  situ 
characterization  of  water  management  related  properties  of  Nation, 
under  a  wide  variety  of  conditions,  i.e.  temperature,  pressure, 
hydration,  and  mechanical  stress,  on  water  transport  across  the 
membrane.  To  obtain  absolute  water  content  imaging,  one  could 
calibrate  image  intensity  of  the  resultant  images  by  simultaneously 
measuring  an  isolated  yet  closely  positioned,  saturated  Nation  layer 
with  known  water  content  as  reference.  In  addition,  the  techniques 
employed  in  this  study  can  be  translated  to  measurement  in  oper¬ 
ational  PEM  fuel  cells. 
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